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FENILFOSFONATO DE CÁLCIO INTERCALADO 
COM SILANO. SÍNTESE, CARACTERIZAÇÃO E 
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Angélica Machi Lazarin 
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RESUMO: O composto fenilfosfonato de 
cálcio (CaPP) foi preparado e intercalado 
como3-aminopropiltrimetoxisilano (3APTS), 
resultando no composto CaPP/3APTS. Os 
compostos foram caracterizados através da 
análise elementar, espectros na região do 
infravermelho, microscopia eletrônica de 
varredura e difração de raios X. Em seguida 
a superfície do composto CaPP/3APTS foi 
adsorvido o fon Cu(Il) (CaPP/3APTS/Cu(Il)). 
A isoterma de adsorção foi ajustada à 
equação modificada de Langmuir, obtendo- 
se a linearização da mesma. Através do 
coeficiente angular calculou-se o valor da 
constante de estabilidade média que foi de 
” = 633 dm* mol” e pelo coeficiente linear 
obteve-se o valor do número médio de 
ligantes imobilizados e simultaneamente 
coordenados a um íon Cu(ll) (7 = 4). O 
valor da capacidade máxima de adsorção 
de Cu(Il) foi de 0,95 mmol.g””. 
PALAVRAS-CHAVE: Intercalação, 
adsorção, 3-aminopropiltrimetoxisilano 
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CALCIUM PHENYLPHOSPHONATE 
INTERCALED WITH SILANE. 
SYNTHESIS, CHARACTERIZATION 
AND ADSORPTION OF CU(II) 

ABSTRACT: The 
phenylphosphonate 
prepared and intercalated with 
3-aminopropyltrimethoxysilane  (3APTS), 
resulting in the compound CaPP/3APTS. 
This compound was characterized through 
elemental analysis, infrared spectra, 
scanning electron microscopy and X-ray 
diffraction. Then, the surface of the 
CaPP/SAPTS compound was adsorbed 
with the Cu(ll) ion (CaPP/3APTS/Cu(Il). 
The adsorption isotherm was adjusted to 
the modified Langmuir equation, obtaining 
its linearization. Through the angular 
coefficient, the value of the average stability 
constant was calculated, which was T = 633 
dm mol and through the linear coefficient, 
the value of the average number of ligands 
immobilized and simultaneously coordinated 
to a Cu(ll) ion (= 4). The value of the 
maximum adsorption capacity of Cu(ll) was 


compound calcium 
(CaPP) was 


0.95 mmol.g”*. 
KEYWORDS:  Intercalation, adsorption, 
3-aminopropyltrimethoxysilane 
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INTRODUÇÃO 


Os compostos inorgânicos cristalinos lamelares se agrupam numa classe que se 
caracteriza por apresentar uma sequência infinita de camadas (ALBERTI, BEIN, 1996). 
Nestes últimos anos, o procedimento experimental com estes tipos de compostos, 
influenciou a Química do estado sólido e a síntese de materiais. Historicamente, o 
impulso inicial deste campo aconteceu, quando pesquisadores voltados à busca de 
trocadores inorgânicos, para serem usados em reatores atômicos, que suportassem altas 
temperaturas, como acontece com os óxidos e seus derivados, obtiveram casualmente 
compostos cristalinos. O primeiro composto cristalino sintetizado foi o fosfato de zircônio, 
que abriu assim, uma corrida ao estudo das propriedades destes materiais, devido ao fato 
de que os mesmos apresentavam espaços vazios entre as lamelas inorgânicas, além do 
fato relacionado à ligação química envolvida na manutenção da estabilidade das camadas 
inorgânicas (ZHANG, CLEARFIELD, 1992). 

Os compostos lamelares são essencialmente de caráter bidimensional no sentido de 
que, as forças de ligações dentro da camada são mais intensas do que aquelas existentes 
entre as mesmas. Assim, uma grande variedade de compostos inorgânicos cristalinos 
possui propriedades estruturais, destacando-se em particular, a formação do esqueleto 
inorgânico em arranjo lamelar. Estes sólidos organizados em camadas desempenham 
papel importante na adsorção, na troca iônica, no funcionamento como peneira molecular, 
na atuação como catalisador, na condução de prótons e íons, além disto, é também 
uma conveniente matriz para modificações químicas. Está associada a estas matrizes 
inorgânicas uma alta estabilidade térmica e relativa resistência química à oxidação e alta 
seletividade a certos íons e moléculas (VIEIRA, LAZARIN, 2017; MIYAMOTO, KURODA, 
OGAWA, 2004). 

As características estruturais destes compostos permitem desenvolver um dos 
mais importantes aspectos inerentes a estes tipos de materiais, que consiste na entrada 
de moléculas neutras e espécies carregadas no espaço definido pela cavidade entre 
as lamelas, provocando uma expansão da distância interlamelar, após a acomodação 
da espécie convidada dentro do hospedeiro inorgânico (CAPKOVÁ, SCHENK, 2008; 
ROCHA, et.al., 2022, LAZARIN, et al., 2022). De um modo geral, estes tipos de matrizes 
são plenamente sensíveis à intercalação, que normalmente ocorre de maneira reversível, 
porém, a estrutura lamelar é mantida com a inserção da espécie convidada. Como resultado, 
é possível manipular compostos lamelares para formar nova e não usuais estruturas. 

Com a inserção de moléculas nos compostos cristalinos lamelares ocorre reação 
interativa com a entrada da espécie no espaço livre da galeria. O efeito contribui para o 
entendimento do comportamento do composto lamelar. Este processo começa com uma 
reorganização do sólido, que ocasiona um aumento na distância interlamelar, que pode ser 
observado através dos picos em baixos valores de 2qg, que é mostrado através da técnica 
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de raios X (MIYAMOTO, KURODA, OGAWA, 2004). 

A caracterização dos compostos intercalados é um grande desafio, sendo 
fundamental para o seu desenvolvimento, não só nos que se refere às propriedades como 
também às aplicações destes sistemas. 

Os fosfonatos intercalados com grupos organofuncionais orgânicos têm despertado 
interesses consideráveis recentemente (RUIZ, AIROLDI, 2003). A organofuncionalização é 
um processo que modifica a superfície de adsorventes porosos, através da ligação química 
com grupos organofuncionais, na qual uma reação entre o suporte e um agente modificador 
adequado, permitindo diversas aplicações e que podem atuar como ligantes para a 
coordenação de íons metálicos na superfície e também na adsorção e pré-concentração 
de íons metálicos de soluções etanólicas (VIERIA, et al., 2017; TRCHOVÁ, et al., 1999). 

Neste trabalho o composto fenilfosfonato de cálcio foi preparado e intercalado com 
o 3-aminopropiltrimetoxisilano. O composto foi caracterizado através da análise elementar, 
espectros na região do infravermelho e difração de raios. Neste composto lamelar foi 
adsorvido o fon Cu(Il), com formação de complexo lamelar metalado (CaPP/3APTS//Cu(Il)). 


MATERIAIS E MÉTODOS 


Síntese do fenilfosfonato de cálcio (CAPP) 


O CaPP foi adquirido misturando-se 60,0 cmº de solução 1,0 mol dm* de ácido 
fenilfosfônico, C,H.PO,H,, com 24,0 cm? de uma solução de CaCl,.2H,0 (VIEIRA, et.al., 
2017). Em seguida, uma solução 1,0 mol dm* de hidróxido de sódio foi adicionada para 
ajustar o pH entre 5 e 6 e para precipitar o sólido branco, o qual foi filtrado e seco a 322 K. 

A reação desta preparação pode ser genericamente escrita como: 

CaCl,.2H,0 + 2[C,H,PO(OH),.., — Ca(HO,PC,H,),.2H,0,, + 2HCl (1) 
Intercalação 

Aproximadamente 45 mg do fenilfosfonato de cálcio em um erlenmeyer com 20,0 
cmº de solução etanólica 1,0 mol dm do 3-aminopropiltrimetoxisilano foi agitado em banho 
termostatizado a 298 + 1 K durante 10 h. O sobrenadante foi desprezado e sólido foi seco 
em estufa a 322 K. 


Caracterização 


As determinações de carbono, hidrogênio e nitrogênio foram feitas em um 
instrumento de Perkin Elmer modelo PE 2400. 

Para a determinação de cálcio e fósforo três amostras de aproximadamente 50 
mg do fenilfosfonato de cálcio foram colocadas em frascos de decomposição de ostefron 
para microondas com capacidade para 90 cm, para digestão ácida, na presença de 2,1 
cmº de ácido nítrico (65 %), 2,1 cmº de ácido fluorídrico (48 %) e 1,0 cm? de peróxido de 
hidrogênio (30 %). Em seguida essas amostras foram levadas a um forno de microondas 
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da Provecto analítica, Modelo DGT 100 plus, para mineralização dos compostos. O tempo 
total utilizado foi de 18 min. Esperou 1 h para a abertura dos frascos de decomposição de 
ostefron e, logo em seguida, o ácido fluorídrico foi eliminado das amostras, na presença 
de aproximadamente 2 g de ácido bórico. Três amostras do composto foram diluídas em 
balão de 250 cmº e posteriormente analisadas. Cada amostra foi medida três vezes através 
do plasma acoplado indutivamente (ICP-AES) da Perkin Elmer 3000DV, por possuir boa 
sensibilidade de detecção multielementar e interferência mínima na amostra analisada. A 
precisão entre as replicatas foi menor que 5 % (AIVES, et al., 2000). 

O aparelho usado na obtenção dos difratogramas de raios-X foi um difratômetro 
marca Shimadzu, modelo XD3A, na faixa de 20 = 2 — 70º, com radiação de Cu Ka. 

A microscopia eletrônica de varredura para os compostos foram efetuadas em um 
aparelho JEOL, modelo JSTM-300 F, de resolução de 1,2 nm em 15 kV e voltagem de 
aceleração de 20 kV. As amostras foram preparadas em um suporte metálico sobre uma 
fita de carbono e metalizadas com uma camada de ouro por 400 s. 

Os espectros de absorção na região do infravermelho foram adquiridos através do 
espectrofotômetro FTIR de marca Bomem, modelo MB-Series, pelo método da pastilha 
auto suportadas em KBr, em 50 varreduras na região entre 4000 e 400 cm", com resolução 
de 4 cm”. 


Isoterma de adsorção 


A adsorção do íon metálico, Cu(ll), na superfície foi efetuada em meio etanólico (25,0 
cm), nas concentrações 2,0 x 102 a 2,0 x 102 mol dm*. Para cada amostra, adicionou-se 
aproximadamente 100 mg de CaPP/3APTS, que foram agitadas mecanicamente por três 
horas, a temperatura ambiente. Após a decantação, pipetou-se 4,00 cmº do sobrenadante 
para a determinação em triplicata das concentrações do metal no equilíbrio (C .), por titulação 
direta com EDTA 0,010 mol dm, usando-se como indicador murexida (ROMANOWSKI, 
MANGRICH, NEVES, 2001; SILVA, FONSECA, ARAKAKI, 2007). Por diferença do número 
de mols do metal adicionado (n,) e o número de diferença do número de mols no equilíbrio 
(n.), dividido pela massa do CaPP (m), determinou-se o valor de (n)), número de mols de 
cobre adsorvido na superfície por grama de CaPP/3APTS, sendo, portanto uma medida 
do grau de cobertura para cada uma das condições experimentais usadas (GURGEL, GIL, 
2009). 

Considerando-se que o íon cobre divalente, quando adsorvido sobre a superfície, 
forma complexos imobilizados com diferentes razões entre ele e o ligante ancorado e 
que a concentração total do ligante imobilizado na superfície é igual ao somatório das 
concentrações de todas as espécies, sendo o número médio de ligantes (n) suportados 
ligados ao íon metálico e a forma linearizada da isoterma é dada pela equação (2) 
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onde T', é o coeficiente de seletividade. 


RESULTADOS E DISCUSSÕES 


A análise elementar de cálcio, fósforo, carbono e hidrogênio para o composto 
sintetizado fenilfosfonato de cálcio foram 10,1, 15,6, 37,2 e 3,0 %, respectivamente. 
Esses valores estão bem próximos daqueles calculados 10,25, 15,90, 36,96 e 3,12 %. A 
partir desses resultados propôs-se a seguinte fórmula molecular Ca(HO,PC,H,),.2H,0. A 
quantidade do 3-aminopropiltrimetoxisilano intercalado no fenilffosfonato de cálcio foi de 
2,3 % de átomos do nitrogênio, que corresponde 1,6 mmol g”. 

Um esquema da intercalação para o 3-aminopropiltrimetoxilano é mostrado na 
Figura 01. 


NH NH, NH, 


NH, NH; NH» 


Figura 01. Representação esquemática do fenilfosfonato de cálcio intercalado com o 
3-aminopropiltrimetoxilano. 


O difratograma de raios X fornece o ângulo entre os raios de incidência e o difratado 
na amostra (20) mostrando ainda a intensidade relativa. A técnica não é destrutiva e apenas 
pequena quantidade da amostra em pó é suficiente para a identificação da fase cristalina 
presente. Através da equação de Bragg, 


2d sen9 = nÃ (3) 


em que | =1,5406 Ão comprimento de onda da radiação utilizada, pode-se determinar 
a distância interlamelar d (À) entre os planos atômicos da rede cristalina. 

O difratograma de raios X do fenilfosfonato de cálcio é mostrado na Figura 02a e 
apresenta picos finos e intensos, caracterizando a cristalinidade do composto. A primeira 
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reflexão aparece em 20 = 5,92º correspondente a uma distância interlamelar referente a 
1494 pm. Como esperado, esta distância aumenta quando o 3-aminopropiltrimetoxisilano é 
intercalado, sendo 20 = 5,52º correspondente a uma distância interlamelar igual a 1601 pm 
como ilustrado na Figura 02b. Um aumento 107 pm em comparação com a matriz original, 
sugerindo que a intercalação ocorre exclusivamente na superfície externa do suporte 


inorgânico (Fig. 01). 


Iiua 


28/erm 


Figura 02. Difratograma do fenilfosfonato de cálcio (a), intercalado com o 3-aminopropiltrimetoxisilano 


(b). 


A Figura 03a apresenta o espectro na região do infravermelho do composto 
fenilfosfonato de cálcio, sendo que a água coordenada ao metal apresenta bandas de 
estiramento e deformação em 3400 cm"! e 1666 cm”!. A banda característica da água é vista 
a 3500 cm”, que se intensifica por causa da fácil adsorção que ocorre durante o processo 
de manipulação da amostra. As bandas em 1439 cm”! e as bandas médias no intervalo de 
721 a 695 cm” são do anel fenílico. As bandas entre 1342 a 1016 cm” são do grupo PO, 
(AIROLDI, ROCA, 1996). As bandas no intervalo de 3000 a 2905 cm são atribuídas ao 
estiramento simétrico e assimétrico de C-H do anel e em 1438 cm” tem-se o estiramento 
C-C do anel (LIMA, AIROLDI, 2001). No espectro do fenilfosfonato de cálcio intercalado 
com o 3-aminopropiltrimetoxisilano (Figuras 03b) não se verifica a presença de bandas 
atribuídas às vibrações do grupo N-H e as outras bandas da matriz inorgânica sofreram 
pequenas mudanças com a intercalação. 
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Figura 03. Espectro na região do infravermelho do CaPP (a) e CaPP/SAPTS (b). 


A Figura 04 mostra a micrografia do fenilfosfonato de cálcio e intercalado com o 
3-aminopropiltrimetoxisilano, respectivamente. Foi observado o aparecimento de placas 
paralelas mostrando que a morfologia do cristal é lamelar (PINTO, PEIXOTO, LAZARIN, 


2022). 


(a) 


Ciências exatas e da terra: teorias e princípios Capítulo 1 


(b) 


Figura 04. Microscopia eletrônica de varredura do fenilfosfonato de cálcio (a) e intercalado com o 
3-aminopropiltrimetoxisilano (b). A barra da escala é dada em micrômetro. 


Isoterma de adsorção 


Uma propriedade do composto cristalino lamelar intercalado é adsorver íons 
metálicos de solução etanólica, portanto, investigou-se a isoterma de adsorção do íon 
metálico cobre (Il). Inicialmente, a investigação demonstrou que a matriz original, sem 
3APTS, não adsorve este íon metálico. 

Na Figura 05 encontram-se a isoterma de adsorção do íon metálico Cu(ll) no 
CaPP/3APTS e na Figura 06 a forma linearizada da isoterma. O número máximo de mols 
adsorvidos (n,) foi 0,95 mmol g! para CaPP/3APTS/Cu(ll). A constante de estabilidade 
média ([,) e o número médio de ligantes ligados e coordenados ao íon metálico (A) foram 
determinados a partir do gráfico de 1/n, e 1//CuCI,. A constante de estabilidade média 
foi 633 dm? mol! e o número médio de ligantes para CaPP/3APTS foi quatro. O valor da 
constante calculada indica a formação de um complexo estável. 
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Figura 05. Isoterma de adsorção do íon metálico Cu?* imobilizado na superfície do CaPP/3APTS, em 
solução etanólica, a 25 ºC. 
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Figura 06. Forma linearizada da isoterma de adsorção do íon metálico Cu?* imobilizado na superfície do 
CaPP/3APTS em solução etanólica. 


CONCLUSÕES 

O fenilfosfonato de cálcio foi preparado e este composto mostrou-se cristalino e 
lamelar comprovados através da difração de raios X e microscopia eletrônica de varredura. 
O resultado das isotermas de adsorção comprovou a obtenção do CaPP/3APTS/Cu(Il). 
A disponibilidade dos centros básicos no composto intercalado confere a este material 
a propriedade de adsorver cátions da soluções etanólicas na interface sólido/líquido, 
comportamento que pode ser útil na remoção de íons metálicos.. 
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RESUMO: O composto fosfato de cálcio 
(CaP) foi preparado e intercalado com o 
ácido p-aminobenzóico (PABA), resultando 
no composto CaP/PABA. Em seguida a 
superfície do composto CaP/PABA foi 
adsorvido o íon Cu(ll) (CaP/PABA/Cu(Il)). 
A isoterma de adsorção foi ajustada à 
equação modificada de Langmuir, obtendo- 
se a linearização da mesma. Através do 
coeficiente angular calculou-se o valor do 
coeficiente de seletividade que foi de T 
= 453 dmº mol! e pelo coeficiente linear 
obteve-se o valor do número médio de 
ligantes imobilizados e simultaneamente 
coordenados a um íon Cu(ll) (à = 4.0 
valor da capacidade máxima de adsorção 
de Cu(ll) foi de 1,75 mmol.g!. Estes 
materiais forão caracterizados por análise 
elementar, difração de raios X, análise 
termogravimétrica e estudo eletroquímico 
utilizando a técnica de voltametria cíclica. 
PALAVRAS-CHAVE: Eletrodo, voltametria 
cíclica, cobre divalente. 
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DEVELOPMENT OF A CARBON 
PASTE ELECTRODE WITH 
DIVALENT COPPER ADSORVED 
AND INTERLACED IN LAMELLAR 
COMPOUND 

ABSTRACT: The 
compound | calcium 
was synthesized and intercalated with 
p-aminobenzoic acid (PABA), resulting 
in the compound CaP/PABA. Then, the 
surface of the CaP/PABA compound was 
adsorbed with the Cu(Il) ion (CaP/PABA/ 
Cu(ll)). The adsorption isotherm was 
adjusted to the modified Langmuir equation, 
obtaining its linearization [ = 453 dm* mol 
t and through the linear coefficient, the 
value of the average number of ligands 
immobilized and simultaneously coordinated 
to a Cu(ll) ion ( = 4). The value of the 
maximum adsorption capacity of Cu(Il) 
was 1,75 mmol g!. These materials will be 
characterized by elemental analysis, X-ray 
diffraction, thermogravimetric analysis and 
electrochemical study using the technique 
of cyclic voltammetry. 
KEYWORDS: Elecirode, cyclic voltammetry, 
divalent copper. 


crystalline lamellar 
phosphate (CaP) 
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INTRODUÇÃO 


O uso de eletrodos quimicamente modificados tem sido amplamente investigada, 
devido às vantagens que esses dispositivos podem fornecer, como a seletividade e 
sensibilidade (EDMONDS, 1988). 

Devido à importância do monitoramento de aplicações bioquímicas e biomédicas, 
o desenvolvimento de novos métodos de análise por meio de sensores eletroquímicos 
baseados em eletrodos de pasta de carbono (KARIMI-MALEH, et al., 2009) fibra de carbono 
(HRBÁC, et AL., 2007), carbono vítreo (EVTUGYN, et al, 2007) e revistidos com polímero 
(Zeng, et al., 2007) são estudados. 

Numerosos compostos inorgânicos são adequados como hospedeiros para este 
tipo de investigação, no entanto, compostos de fosfato cristalinos lamelares apresentam 
alta capacidade de troca iônica e alta condutividade (CLEARFIELD, 1988), características 
que são muito atraentes para utilizar este material como uma matriz para imobilização de 
espécies eletroativas (RAJ, OHSAKA, 2001; LAZARIN, et al., 2022). 

Neste trabalho o composto fosfato de cálcio foi preparado e intercalado com o 
ácido p-aminobenzóico. Em seguida a superfície do composto CaP/PABA foi adsorvido o 
fon Cu(Il) (CaP/PABA/Cu(Il)). Os compostos forão caracterizados por análise elementar, 
difração de raios X, análise termogravimétrica e estudo eletroquímico utilizando a técnica 


de voltametria cíclica. 


MATERIAIS E MÉTODOS 


Síntese do fosfato de cálcio (CAP) 

O fosfato de cálcio foi preparado adicionando-se vagarosamente uma solução 0,50 
mol dm?* de CaClI,.2H,0 a uma solução 1,50 mol dm-? de fosfato de amônio dibásico e 
aquecida à 363 K (SILVA, et al., 2012). A suspensão formada permaneceu reagindo durante 
1 he em seguida o precipitado obtido foi filtrado e seco à 324 K. Em seguida o composto foi 
aquecido à 433 K, durante 48 h, para eliminar a amônia, as reações foram escritas: 


CaCI1,.2H,0 2(NH,),HPO, — (NH,),Ca(HPO,),.H,0 + 2NH,CI+H,O 
(NH,),Ca(HPO,),.H,0 — Ca(H,PO,), + 2NH, + H,O 


Am 
1º 
0 — 


Intercalação 


Em torno de 3,0 g do fosfato de cálcio num erlenmeyer com 25,0 cm* de solução 
etanólica 1,2 mol dm do ácido p-aminobenzóico foi agitado em banho termostatizado a 
298 + 1 K por 12 h. O sobrenadante foi desprezado e sólido foi seco em estufa a 323 K. 
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Caracterização 


As análises elementares de carbono, hidrogênio e nitrogênio foram feitas num 
instrumento de Perkin Elmer modelo PE 2400. 

As análises elementares de cálcio e fósforo (Alves, et al., 2000) foram realizadas 
por espectroscopia de absorção atômica utilizando um espectrômetro de absorção atômica 
Perkin Elmer, modelo 5100, e métodos espectrofotométricos utilizando um espectrofotômetro 
Shimadzu, modelo MultiSpec-1501. 

O aparelho utilizado na obtenção dos difratogramas de raios-X foi um difratômetro 
marca Shimadzu, modelo XD3A, na faixa de 20 = 2 — 70º, com radiação de Cu Ka. 

As curvas termogravimétricas foram obtidas em um instrumento da DuPont, modelo 
1090 B, acoplado a uma termobalança 951, em intervalo de temperatura compreendido 
entre 273 a 1273 K, com taxa com razão de aquecimento de 0,167 K s?, sob atmosfera 
dinâmica de argônio de 1,66 cm s'!, com massa de aproximadamente 10 mg em todas as 


amostras. 
Isoterma de adsorção 


Para obtenção da isoterma de adsorção de cobre pela CaP/PABA, preparou-se 
inicialmente uma solução estoque de CuCIl, 0,1 mol dm* em etanol anidro. Diferentes 
volumes dessa solução foram pipetados e diluídos a 25,0 cmº com etanol anidro, a fim de 
se obter amostras com concentrações definidas (C) que variaram de 1,0 x 10% a 1,0 x 102 
mol dm. 

Acada amostra adicionou-se cerca de 0,100 g do composto intercalado e estas foram 
agitadas por quatro horas, a temperatura de 298 K. Pipetou-se 6,0 cmº do sobrenadante, 
para determinar as concentrações do metal no equilíbrio (C ), realizadas por titulação direta 
com EDTA (0,010 mol dm), usando-se murexida como indicador (BENETTI, et al., 2022). 

Pela diferença do número de mols de Cu(Il) adicionado (n.) e o número de mols de 
metal no equilíbrio (n,), dividido pela massa do CaP/PABA (m), determinou-se os valores 
de (n,), número de mols de metal adsorvido na superfície por grama do CAP/PABA, sendo a 
medida do grau de cobertura para cada uma das condições experimentais utilizadas, isto é: 


Cu(IN)] 


m 


(3) 


h, = 


m 


mom [(CaP/ PABA) 
m m 


visto que a somatória estende em todas as espécies sobre a superfície. O número 
máximo de mols do complexo adsorvido n”* é igual ao valor de n, quando a concentração 
de Cu(Il) tende a infinito. 

A adsorção do íon cobre divalente contido em solução, depende dos coeficientes de 
seletividade dos complexos formados com o ligante imobilizado na superfície. O equilíbrio 


de formação do complexo de cobre divalente com o composto intercalado foi estudado com 
base no modelo de Filippov (FILIPPOV, 1984; ROCHA, et. al, 2022). 
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A reação de formação do complexo na superfície pode ser representada 
genericamente pela equação 4: 


MY,+nR o MRY, (4) 


sendo que z é a carga do íon metálico, Y é um ânion com carga unitária e Ré o 
grupo ligante, eletricamente neutro, imobilizado na superfície do suporte. 
No equilíbrio (5), pode ser aplicado uma relação do tipo: 


9,/(1-0,)= y,C (5) 


sendo que 6, é a fração de ligante ligado ao metal, sendo que 6 =n n/N;Néa 

concentração total do ligante imobilizado na superfície (mol/g); n é o número de ligantes 

suportados no composto intercalado e coordenado ao íon metálico; y, é o coeficiente de 
seletividade do complexo suportado. 

Lo à ê Ã 

ny, Ny, N 


(6) 


Considerando-se que o íon cobre divalente, quando adsorvido na superfície, pode 
formar complexos imobilizados com diferentes razões entre ele e o composto intercalado 
ancorado, que a concentração total do ligante imobilizado na superfície é igual ao somatório 
das concentrações de todas as espécies; (/i) o número médio de ligantes suportados 
ligados ao íon metálico será: 


à = UMRY JH 2IMRSY. +... m[MR,Y.] 


m z (7) 
e 
A equação 6 pode ser escrita: 
a = n + Ea (8) 
n, MSC N 


onde T, é o coeficiente de seletividade iguala > Y,, 


m=1 

Estudo eletroquímico 

Medidas eletroquímicas foram realizadas em um Potenciostato 263A.EG&C PAR. O 
eletrodo de referência usado foi o eletrodo Ag/AgCI saturado com KCI, dupla junção, Orion, 
modelos 90-02 (Eo = 0,197 V vs EPH) e o eletrodo auxiliar foi um fio de platina. O eletrodo 
de trabalho utilizado para o estudo foi o eletrodo de pasta de carbono CaP/PABA e CaP/ 
PABA/Cu(lIl). As pastas foram feitas misturando-se 25 mg do composto, 35 mg de grafite e 
uma gota de nujol. 

Os eletrodos foram confeccionados colocando-se na extremidade de um tubo de 
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vidro de 5 mm de diâmetro interno uma placa de platina, deixando uma cavidade de 2 mm 
na extremidade para colocação da pasta. A esta placa foi soldado um fio de platina e, a 
este, um fio de cobre para contato elétrico com o aparelho. 

Os eletrodos foram confeccionados colocando-se na extremidade de um tubo de 
vidro de 5 mm de diâmetro interno uma placa de platina, deixando uma cavidade de 2 mm 
na extremidade para colocação da pasta. A esta placa foi soldado um fio de platina e, a 
este, um fio de cobre para contato elétrico com o aparelho (Figura 01). 


fio de Cobre 


tubo de vidro 


fio de Platina 


placa de Platina 
pasta de Grafite 


Figura 01. Esquema do eletrodo de trabalho. 


RESULTADOS E DISCUSSÕES 


As quantidades de cálcio e fosfóro para o fosfato de cálcio foram 26,4 e 17,0 %, 
respectivamente. Esses valores estão bem próximos daqueles calculados 26,5 e 17,1 %. A 
fórmula molecular Ca(H PO), foi proposta a partir dos resultados encontrados. 

A quantidade de PABA intercalado no CaP foi de 6,35 % de átomos do nitrogênio, 
que corresponde 4,54 mmol g!. O mecanismo de intercalação do PABA nas camadas 
inorgânicas pressupõe a protonação dos centros básicos pelos grupos ácidos P-OH 
localizados na matriz inorgânica, em uma típica interação ácido-base.O esquema da 
intercalação para o PABA dentro da cavidade inorgânica é mostrado na Figura 02. 
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Figura 02. Representação esquemática do fosfato de cálcio intercalado com o ácido p-aminobenzóico. 


Os difratogramas de raios X do fosfato de cálcio e intercalado com o ácido 
p-aminobenzóico apresentam picos finos indicando uma boa cristalinidade. O pico em 
20 = 12,76º relativo à primeira reflexão referente ao plano (002), informa que à distância 
interlamelar do fosfato de cálcio é igual a 697 pm (AIROLDI, LIMA, 2004). A distância 
interlamelar aumenta para 1570 pm quando o fosfato de cálcio é intercalado com o ácido 
p-aminobenzóico. Um aumento de 873 pm em comparação com a matriz original, mostrando 
que a intercalação ocorre na superfície interna do suporte inorgânico (Fig. 02). 

As curvas termogravimétricas do CaP e intercalado com o PABA apresentam 
perda de massa em 635 K que corresponde a perda da amina quimicamente adsorvida e 
duas moléculas de água para a formação de pirofosfato de cálcio (MURAKAMI, 1992). A 
curva termogravimétrica para o CaP/PABA é idêntica ao do CaP, mostrando diferença na 
percentagem no estágio de decomposição, sendo maiores para o composto intercalado. 


Isoterma de adsorção 


Na Figura 03 encontra-se a isoterma de adsorção do fon metálico e na Figura 04 a 
forma linearizada da isoterma. O valor obtido para o número máximo de mols adsorvidos 
(n) foi de 1,75 mmol g”', o número médio de ligantes coordenados ao íon metálico (1 ) 4 e 
o coeficiente de seletividade ([,) 453 dmº.mol”. 
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Figura 083. Isoterma de adsorção do íon metálico Cu?* imobilizado na superfície do CaP/PABA, em 
solução etanólica, a 25 ºC. 
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Figura 04. Forma linearizada da isoterma de adsorção do íon metálico Cu?* imobilizado na superfície do 
CaP/PABA em solução etanólica. 


Estudo eletroquímico 

Após a preparação dos eletrodos quimicamente modificados foram obtidos os 
voltamogramas cíclicos, no intervalo de potencial de -1,0 a 1,0 V, sob argônio. Para o 
eletrodo modificado com CaP/PABA nenhuma corrente de pico foi observada, entretanto, 
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um acoplamento redox foi verificado ao se usar o eletrodo CaP/PABA/Cu(Il). Nesse caso, 
observou-se que o potencial médio foi de E, = 0,16 V, sendo E, = (Ep, + Ep.)/2; Ep, e 
Ep, são os potenciais de pico catódico e anódico, respectivamente. Este potencial de pico 
anódico perto de zero é muito importante para o desenvolvimento de sensor/biosensor. 
Em comparação com alguns resultados obtidos, esta onda pode ser atribuída ao seguinte 
processo de transferência de elétron (BORGO, et al., 1999): 


[Cu"(CaP/ABA),|P* + e —— [Cu(CaP/ABA) |) (9) 


A verificação da estabilidade química do complexo de cobre divalente sobre o CaP 
intercalado com o PABA, foram feitas várias varreduras numa faixa de potencial entre -1,0 
e 1,0 V vs ECS a uma velocidade de varredura de 20 mVs'! com o eletrodo modificado 
imerso numa solução de tampão fosfato 0,11 mol dm. Os resultados mostram que não 
há variação significativa antes de 100 ciclos redoxes, como mostrados na Figura 05. As 
intensidades de corrente de pico permaneceram constantes, logo o complexo de cobre 
divalente está aderido sobre o CaP intercalado com o PABA. 
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Figura 05. Gráfico de corrente de pico, Ip vs número de ciclos, para o eletrodo CaP/ABA/Cu(Il), sob 
atmosfera de argônio.Eletrólito suporte: tampão fosfato 0,10 mol dm*?. 


A resposta de corrente de pico do eletrodo CaP/ABA/Cu(ll) em vários ciclos 
voltamétricos apresentou uma dependência em relação ao aumento da velocidade de 
varredura. Através do gráfico Ip vs. v!?, observou se uma dependência linear (Figura 06), 


que é característica de um comportamento similar àqueles associados com o processo 
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difusional (PANICE, et al., 2014). Sistemas eletroquímicos controlados por difusão, com 
cinética da reação limitada pela difusão das espécies do interior da solução até a superfície 
do eletrodo (difusão semi-infinita) apresentam uma relação linear das correntes de pico (lp) 
em função a raiz quadrada da velocidade de varredura (lp x v!2) (KUBOTA, et al., 1995). 

A natureza do eletrólito suporte não produziu nenhuma influência significativa no 
E, cujos valores variaram entre 0,16 e 0,17 V. Os resultados mostraram que o eletrólito 
suporte interage com a superfície da matriz e não produz nenhuma mudança no potencial 
médio. 


DAD D.15 D.2D DZ5 D,.ID 


Figura 06. Gráfico de corrente de pico, Ip vs v!?, para o eletrodo CaP/PABA/Cu(Il), sob atmosfera de 
argônio. Eletrólito suporte: tampão fosfato 0,10 mol dm”. 


A oxidação do ácido ascórbico na superfície do eletrodo CaP/ABA/Cu(Il) é mostrado 
na Figura 07. A curva de voltametria cíclica na presença de 1,0x10% mol dm* de solução de 
ácido ascórbico em pH 7,2 mostrou um considerável aumento da corrente de pico anódico 
do eletrodo. 
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Figura 07. Voltamograma cíclico do eletrodo CaP/PABA/Cu([Il) na ausência (a) e na presença de ácido 
ascórbico (b) (solução tampão fosfato 0,10 mol dm) pH 7,2. Velocidade de varredura: 20 mVs”. 


CONCLUSÕES 


O CaPfoiobtido e este apresentou-se cristalino e lamelar. Através da termogravimetria 
observaram-se perdas de massas relativas à saída de água e deixando o resíduo de 
pirofosfato de cálcio e no fosfato de cálcio intercalado com o ácido p-aminobenzóico. 
Quando o CaP é intercalado com o PABA há elevação da distância interlamelar para 
acomodação das moléculas 

A análise do resultado da isoterma de adsorção forneceu subsídio comprobatório da 
obtenção da CaP/PABA adsorvida com cobre divalente. O complexo de cobre divalente foi 
imobilizado no CaP/PABA. Não observou lixiviação do complexo.. A natureza do eletrólito 
suporte não afetou significativamente a resposta do eletrodo. O potencial de pico anódico 
sob o qual o ácido ascórbico é oxidado utilizando-se o eletrodo CaP/PABA/Cu(Il) é baixo. 
O eletrodo não apresentou mudanças em suas respostas após 100 determinações, o 
que caracteriza uma boa estabilidade quando se trata de eletrodo de pasta de carbono. 
O eletrodo é estável sob as condições de operação. Os resultados obtidos usando o 
presente eletrodo são muito bons considerando que outras espécies químicas interferem 
nos resultados. As vantagens em usar este eletrodo são principalmente a facilidade com 
que ele é preparado e também o pouco tempo de resposta que é gasto no procedimento 
analítico. Essas características fazem com que o presente material seja muito atrativo para 


ser utilizado como sensor para ácido ascórbico. 
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ABSTRACT: The porosity and apparent 
density of iron ore pellets directly interfere 
with the blast furnace process and, therefore, 
need to be known to assist in its control 
and optimization. These characteristics are 
generally calculated using a pycnometer 
that uses mercury under pressure to fill 
the pores of the pellet. Considering the 
need to preserve the environment and the 
safety of operators, proposals were made to 
replace this process, but there are several 
complaints about the repeatability of results 
achieved, in addition to the time spent in 
preparing and executing these essays. At 
the same time, it is possible to observe a 
remarkable development in Computer Vision 
and Artificial Intelligence, mainly through 
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Convolutional Neural Networks, which 
can extract patterns from a set of images 
and detect these same patterns in images 
subsequently exposed to this network. In 
addition to performing classification and 
detection, the Mask R-CNN network can 
perform pixel-by-pixel segmentation of 
objects in images. In its evaluation, the 
network presented a significantly high mAP 
and accuracy, demonstrating a satisfactory 
result for the segmentation and obtaining of 
porosity and apparent density values, with 
results similar to the essays currently used. 
KEYWORDS: Convolutional Neural 
Networks, Apparent Porosity, Computer 
Vision, Mask R-CNN, Artificial Intelligence. 


11 INTRODUCTION 


Technological innovations make 


industrial processes increasingly agile, 
fast, efficient, and sustainable. In this 
context, the conservation of efficiency 
and sustainability in processes over time 
becomes crucial for companies. 

The creation of ISO 14000, which 
deals with environmental management and 


auditing, and its required implementation in 
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Brazilian iron ore mining companies by customers of this product by the main foreign steel 
companies, is a reflection that environmental management is increasingly important, not 
only for the important environmental side but also for the economic one as stated by Oliveira 
et al. [1]. 

The porosity and apparent density values of iron ore pellets are key factors in the 
agglomeration and burning processes, belonging to the production stages of iron ore pellets, 
directly impacting the quality of these pellets. This same importance is analyzed when the 
pellets are used as raw material in a blast fumnace process, directly interfering in it, as 
quoted by Fonseca [2]. 

Bayao et al. [3], described in the publication of their study on the influence of porosity 
on the quality of iron” ore pellets, that this characteristic is a parameter that must be well 
controlled during the process, as this property has a strong influence on the index of 
reducibility and compressive strength, which are critical factors for the quality of pellets. 

Still on the importance of the degree of porosity and its relationship with strength, 
Bayao et al. [3] claim that” the strength of an iron ore pellet decreases with increasing 
porosity and if this porosity does not meet the criteria of the steel mills, the pellets can 
break during transport and loading on the blast furnace, reducing its permeability during the 
process and influencing the quality of the steelworks' final product. 

The use of mercury in the process of checking the apparent porosity of an iron ore 
pellet, using a mercuryfilled pycnometer, as described in ISO 15968 and JIS M 8719, does 
not follow current global needs for the environment, making the handling, reuse, packaging 
and disposal of pellets that have gone through these tests a complicated, dangerous and 
expensive task for companies that practice it in this way, as Y. Omori [4] mentions. 

Considering the need to preserve the environment and even due to the development 
of new technologies, this process was replaced by others that do not use mercury, but other 
chemical elements and also water, but there are several reports and complaints regarding 
repeatability of results achieved in these tests using different elements, due to these being 
able to achieve different pore depths in the pellets as stated by Silva [5]. 

Ineffective results are also achieved when an operator, with the aid of a caliper, 
calculates the average diameter of the pellet with ten direct measurements and obtains 
its weight through a precision scale, reaching through these values and others previously 
known, the desired apparent density of the iron ore pellet, as also said by Silva [5]. 

While we observe the evolution of the tests used to analyze the porosity and apparent 
density of iron ore pellets, it is possible to observe a remarkable development in the field of 
Computer Vision and Artificial Intelligence, mainly through Convolutional Neural Networks. 
Through previous training, these Neural Networks can extract patterns and shapes from 
a set of training images and subsequently perform the detection of these same patterns, 
in images exposed to this network, classifying, and detecting objects present in them, 
according to Gonzaga and de Almeida [6]. 
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The R-CNN Neural Mask Network developed by He et al. [7], increments the Faster 
R-CNN network of Ren et al. [8], developing an object mask in parallel to a bounding box, 
performing pixel-by-pixel segmentation of objects in images. The created mask fills the 
pixels of each detected object, creating segregation of it, with the others that are in the same 
static image or analyzed video frame. 

The work proposes to present the use of a previously trained Mask R-CNN network 
to perform the segmentation of iron ore pellets. Through this segmentation and the weight 
of each pellet, achieve the porosity and apparent density values, with results similar to 
the tests described above, bringing more reliability to the processes of iron ore mining 
companies and their client companies, the steel mills. 

This article is divided into five chapters. The first chapter contains a summary of the 
contextualization of the proposed problem. The second presents the theoretical framework 
used to support this research. Chapter three describes the methodology used in this 
research, describing the activities developed and the parameters used for its subsequent 
reproduction. In the fourth chapter, the results achieved will be presented and in the last 
chapter, the conclusions of the work will be offered, with the evaluation of the results and 
proposals for future work. 


21 THEORETICAL REFERENCES 


2.1 Porosity and bulk density of iron ore pellets 


Iron ore is the raw material used in the pig iron production process in a blast furnace. 
This iron ore is supplied in three different forms to steel companies: in the form of Sinter 
Feed, which goes through the sintering process and then goes to the blast furnace, the AF 
Granulate that goes straight to the process, and finally, the Pellet Feed, which goes through 
the entire pelletizing process, to be fed into the blast furnaces of the steel mills. 

During the pelletizing process, the Pellet Feed and other agglomerating materials 
and mixed additives will go through several stages until reaching the pelletizing disc, in 
which they will take the form of small spherical pellets, which vary in size according to what 
is requested by the client steel companies, generally, around 8 to 18mm in diameter Valer 
[9]. 

After the pelletizing stage, the formed pellets are sent to the oven to be baked, 
removing much of the moisture gathered in the previous steps. During pellet cooking, the 
spaces previously occupied by moisture and air give way to “empty” spaces called pores. 

The pores of iron ore pellets are “empty” or gas-filled spaces that may or may not be 
in contact with the surface of the material. When these spaces are not in contact with the 
surface, they are called closed pores and when they are, they are called open pores. Open 
pores are the most harmful to mechanical strength, as many fractures start from superficial 
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defects, as demonstrated by Bayao et al. [3]. Due to the decrease in pellet strength” with the 
increase in the number of pores and its importance for transport processes and even for the 
blast furnace process, the porosity value of the pellets must be determined. The percentage 
occupied by pores in the volume of a pellet can be calculated using eg. (1), informed by 
Silva [5]. 


Apparent Densit 
Pe ioxdi=l pparent Density 


Real Density ) (1) 


As well as the porosity, the Apparent Density of iron ore pellets also directly influences 
the process of a blast furnace and likewise needs to be checked before being supplied. 
According to Silva [5], the Apparent Density of an iron ore pellet is governed by eg. (2). 
Weight 
Volume. (2) 
The same author also adds that the Real Density corresponds to the real volume 


Apparent Density = 


occupied by this pellet, disregarding the value of its porosity. The value of this density is 
based on the iron content of the batch of minerals and binders that form the pellet. 


2.2 Convolutional Neural Networks and Their Further Developments 


The Convolutional Networks of LeCun et al. [10], also known as Convolutional Neural 
Networks, or CNN, is a specialized type of Deep Neural Network, for processing data that 
can be formatted through unidimensional or multidimensional arrays such as images, which 
can be understood as a grid of pixels in two dimensions. In LeCun et al. [10], a multilayer 
Artificial Neural Network was developed, called LeNet-5 that could classify handwritten 
digits from small representations of patterns found in the images presented to the network. 

Several developments after the work of LeCun et al. [10] are presented by Gu et 
al. [11]. The author mentions that since 2006, several methods have been developed to 
overcome the difficulties encountered in training deep CNN, but he mentions that the AlexNet 
Krizhevsky and Sutskever [12] network was the one that triggered the first improvements. 
This network was similar to LeCun et al. [10] LeNet-5 network, but with more layers. Also, 
according to Gu et al. [11], the ResNet network from He et al. [13], that won the ILSVRC 2015 
competition is 20 times deeper than the AlexNet network from Krizhevsky and Sutskever 
[12] and 8 times deeper than the VGGNet developed by Simonyan and Zisserman [14]. The 
addition of layers to CNNs can improve the approximation of the network activation function 
to the target with greater non-linearity, achieving better representative characteristics of 
objects in the images presented to the network. 

Girshick et al. [15] proposed a new CNN configuration in order to extract different 
regions from the same image. This configuration (Fig. 1) was named as Region-based 
Convolutional Neural Network (R-CNN), as it extracts from the analyzed images, around 
2000 regions of interest, instead of trying to classify a huge amount of regions. Subsequently, 
these regions are analyzed by a CNN and the classification of each one of them will be 
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performed using supervised learning models (SVM). With the new configuration developed, 
the authors achieved results of 53.3% of mean Average Precision (mAP), a result that was 
more than 30% higher than the best value obtained by the other competitors of the VOC 
2012. 


 acroplane? no. 
> Eno? yes. 
n E 
4 4 tvmonitor? no. 
1. Input 2. Extract region 3. Compute 4. Classify 
image  proposals (=2k) CNN features regions 


Figure 1. R-CNN model steps. Source: Girshick et al. [15]. 


The next work from Girshick [16] uses his previously mentioned work and others to 
obtain greater efficiency in object detection using Deep Convolutional Neural Networks. 
The previous model has several disadvantages, such as training performed in several 
stages, training with high processing cost, and high memory usage, in addition to slow 
object detection. The related disadvantages were solved by replacing the ConvNet that 
was renewed for each object analyzed by a SPPnet (Spatial Pyramid Pooling networks) 
that used a model of shared characteristics between the analyzed objects. Concerning his 
previous work, the new model from Girshick [16] employed several innovations to improve 
the training and testing speed, in addition to the accuracy in detecting the analyzed objects. 
This new model (Fig. 2) performed the training of the VGG16 network nine times faster and 
achieved higher mAP values in the VOC 2012 challenge, than his previous work. 


Rol feature 
feature map VECLOr  roreachro 


Figure 2. Fast R-CNN model architecture. Source: Girshick [16]. 


Ren et al. [8] made improvements to the Fast R-CNN model from Girshick [16] in order 
to achieve realtime object detection, thus developing the Faster R-CNN model. A large part 
of the reduction in the processing cost was achieved through the sharing of characteristics 
between the networks used and that was part of the Fast R-CNN set. This model almost 
achieved real-time detection, when it ignored the time spent creating the regions that would 
be analyzed, making this the bottleneck of detection systems. To deal with this bottleneck, 
the authors developed a network called Region Proposal Network (RPN) which intended 
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to share the convolutional layers with the networks that performed object detection. During 
the development of the work, the authors observed that the feature maps used by region 
detectors, such as the Fast R-CNN, could also be used to propose regions. Thus, with the 
inclusion of convolutional layers, the edges of each region were reduced, achieving a better 
result in object detection. 

Taking advantage of the architecture of the Faster-RCNN model, He et al. [7] 
developed a new model capable of performing the isolation of objects from images 
presented to the network through filled pixel to pixel masks. This new model was known as 
Mask R-CNN and was used to develop this work. 


2.3 Mask R-CNN 


The method developed by He et al. [7], called Mask R-CNN, increments the Faster 
R-CNN network from Ren et al. [8] by developing an object mask in parallel with a bounding 
box. The mask created fills the pixels of the detected object, creating segregation between 
the others that are in the same static image or analyzed video frame. Although the Faster 
R-CNN by Ren et al. [8] is not designed to have a pixel-level alignment between the inputs 
and outputs of the network, a simple and quantization-free layer, called a RolAlign, which 
faithfully preserves exact spatial locations. Although simple, this change can directly impact 
the improvement of mask accuracy with values between 10 to 50%. With allthe improvements 
implemented, the new method outperformed all competitors in the 2016 COCO challenge, 
reaching rates of 5fps when running on a GPU (200ms per frame of video). 

According to He et al. [7], the object detection system in images is formed by four 
distinct modules. The first of these modules perform the extraction of characteristics of 
each object in the images. The second is a CNN that creates Regions of Interest based 
on the outputs of the previous module, that is, through the features that were extracted 
previously. In the third module, the objects present in each Region of Interest generated in 
the second module are classified by a CNN. The last module is responsible for performing 
the segmentation of image objects, through a binary mask applied to each object. 


31 METHODOLOGY 


For the application of deep convolutional neural networks for the apparent density 
analysis of iron ore pellets, a network based on the Mask R-CNN architecture by He et al. [7] 
was developed, using the learning transfer technique, and fine-tuning from Microsoft COCO 
dataset weights. The weights from this dataset were used as a starting point for developing 
our network weights during training, due to their proven robustness and popularity. 

No dataset was publicly found for iron ore pellets, so it was necessary to develop a 
proprietary dataset containing 44 image files of iron ore pellets, presented individually or 
in groups. The image folder for training the neural network contains 36 images of pellets 
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taken from the internet. The folder used to validate the results, on the other hand, contains 
8 images different ffom those described above, 5 of which belong to the samples that were 
used for comparison with the results of equipment in use, responsible for carrying out porosity 
and bulk density tests of a metallurgical laboratory of a mining company. The resolution of 
images ranges from 259x194 to 4000x2000 pixels. On each folder, an annotation file in 
JSON format was included, containing information on the positioning and segmentation of 
the iron ore pellets contained in each image used for training and validation. The annotations 
were performed using version 2.0.10 of the online tool, VGG Image Annotator (VIA) Dutta 
etal. [17]. 

For the development of the neural network algorithm, Google Collaboratory (Colab), 
an online platform based on the Jupyter Notebook tool, was used. This choice was due to 
the ease in implementing packages developed by third parties and access to the files of 
the neural network created, which could be cloned directly from GitHub servers and were 
hosted on Google Drive. In addition to these facilities, the tool provides limited access to 
segregated Google server GPUSs for this purpose. At the time this work was written, most 
servers used Nvidia Tesla T4 model GPUs, which mainly allowed less time for training and 
making the necessary inferences during this development. 

For the training of the developed Neural Network, 2000 epochs were performed with 
the model, which contains 181 instances of pellets in its dataset. Due to the small number 
of images in the originally created dataset, data augmentation techniques were used to 
increase the number of samples for network training. The techniques used subsequently 
allowed an increase in accuracy from 95% to values close to 98%. 

The pellet samples (5 units) and the porosity reports containing the results of the 
tests carried out on each of these were provided by the Vale metallurgical laboratory, located 
in the city of Vitoria/ES, and the results were” obtained using the equipment called manual 
porosimeter, manufactured by ControlVix. Sample data are shown in Table 1. 


Sample Diameter Diameter Volume Apparent Density Porosity 
(pixels) (mm) (cm?) (g/em?) (%o) 
1 2,425 11,28 0,7508 3,2297 35,82 
2 2,434 11,31 0,7582 3,2103 36,21 
3 2,911 12,08 0,9224 3,1557 37,29 
4 2,181 11,80 0,8604 3,2322 35,77 
5 3,121 12,23 0,9585 3,2562 35,29 


Table 1. Sample data provided. 


After the detection and segmentation of the pellets present in each image exposed 
to the network (Fig. 3), diameter measurements were performed in pixels for each pellet. 
This value was transformed into millimeters, to obtain the actual diameter of the samples. 
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Figure 3. Pellet samples segmented by the Mask R-CNN network. 


To obtain the value in millimeters represented by each pixel of the image, the image 
of a pellet with exactly 12mm in diameter was printed on a blank sheet. This printed image 
was presented to the network at a distance of 10cm from the camera lens, this being the 
same distance for obtaining the images of the samples. The image pellet was detected by 
the net, presenting a mean result of 0.117887mm/pixel. 

Thus, this value was used to perform the calculations of apparent density and then 
porosity, through the formulas previously mentioned, so that it would be possible to directly 
compare the results with the equipment in use in the laboratories of the mining company in 
question. The weight values of each sample provided were checked using a precision scale 
and all these were confirmed. 

To verify the results of the model applied to the study, the mean Average Precision 
(mAP) values will be used, defined by Tan [18], in addition to the accuracy and loss values. 
These values were automatically calculated by the proposed model during its execution. 


41 RESULTS 


The execution of the 2000 epochs of training the model with the dataset described 
in the previous topic, were executed in 15.5 minutes on the Google Colab GPUs and the 
results of accuracy, loss, false positives and false negatives achieved, after the execution of 
the training, were quite favorable for the developed configuration (Fig. 4). 

During the validation of images containing more than one pellet, the algorithm found 
mainly the pellets that were in the foreground and with their outline completely presented 
in the images. At first sight, this non-detection may represent an obstacle, but for this work, 
as we want to present and detect only one pellet in each image presented to the developed 
network, we will not have this concern. 

As seen in the validation image example (Fig. 5), the pellets were segmented 
correctly and their diameter in pixels was calculated and presented by the network. 

With the detection of the pixel diameter results obtained for each of the iron ore pellet 
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samples presented to the network, it was possible to calculate the actual diameter of each 
pellet in millimeters, then eg. 1 and 2 were 


0,9766 | 


% of total 


epochs 


(a) [SR | to) 


Figure 4. Graphs of (a) accuracy, (b) False Positives and (d) False Negatives of the network. 


al " 


102.471 


Figure 5. Image of sample segmented by the Mask R-CNN network, with diameter information in pixels 
of an iron ore pellet sample. 


applied to obtain the percentage of Porosity and Apparent Density of the analyzed 
samples. The actual density value used to calculate Porosity was informed by the supplier 
of the samples, corresponding to 5.03226 g/cm. The volume value was calculated using the 
sphere volume equation. The results achieved were presented in Table 2 and the relative 
errors in Table 3. 


Sample Diameter Diameter Volume Apparent Density Porosity 
(pixels) (mm) (cm?) (g/em?) (%o) 
1 95.455 11.25 0.7461 3.2503 35.41 
2 95.739 11.29 0.7528 3.2334 35.75 
3 102.471 12.08 0.9230 3.1539 37.33 
4 99.946 11.78 0.8564 3.2472 35.47 
5 103.543 12.21 0.9523 3.2774 34.87 


Table 2. Results Achieved Using the R-CNN Mask Network 
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Sample Diameter Volume Apparent Density Porosity 


(mm) (cmê) (g/em?) (%) 
1 -0.24% -0.63% 0.64% -1.14% 
2 -0.21% -0.71% 0.72% -1.28% 
3 0.00% 0.06% -0.06% 0.10% 
4 -0.15% -0.46% 0.46% -0.83% 
5 -0.19% -0.65% 0.65% -1.18% 


Table 3. Relative Errors 


5 | CONCLUSIONS 


Even with a small number of images, the developed dataset, allowed the Mask 
R-CNN network model to be trained to detect, segment, and correctly calculate the diameter 
of iron ore pellets that were exposed to the network, through static images, present in the 
dataset of validation. The network met the expectations, segmenting all pellet samples with 
more than 99% of accuracy during the validation. All relative errors measured were below 
the 1.3% range (Table 3). 

For future work, it is suggested to expand the training and validation dataset, with 
the use of very different samples in format and size. For this expansion, images with 
different lighting and presentation should also be used, increasing the number of unique 
characteristics of the pellets, which could be learned by the model during training. 

tis believed that the error presented in the pellet diameter values (Table 3) comes 
from the difference in methodology as the tests are performed on the equipment in use at 
the mining company. This equipment performs a simple average of the pellet diameters, 
through 60 images of each sample, distributed in 6 different positions. For future work, itis 
suggested that the calculation of the pellet diameter should also be an average of several 
positions of each sample since the pellets do not have a spherical shape that guarantees 
the calculation through only one image. 
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ABSTRACT: Fires are causing major 
disturbances to forests in the Amazon. Its 
cause is related to the interaction between 
the climate and changes in land use and 
cover. A useful tool in this process is the 
quantification of spatial data, which through 
mapping can provide information for 
planning prevention policies. For this, data 
on land use and cover and deforestation 
were used from the MapBiomas project 
for APA Triunfo do Xingu between 1985 to 
2019. The meteorological variables and fires 
were obtained from the Instituto Nacional 
de Meteorologia and BDQueimadas. With 
the results, the loss of vegetation cover 
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over 30 years was approximately 4030.73 
km2 (23.96% of the territory) and for 2004- 
2007 6.21% (1044.84 km2), 2008-2011 of 
3.62% (608.52 km2), 2012-2015 of 1.40% 
(236.24 km2). We can conclude that the 
fires are directly related to the replacement 
of vegetation cover by pastures in the APA 
Triunfo do Xingu, generating a great impact 
on biomass losses. Changes in land use 
were also responsible for a change in air 
temperature, increasing over the past few 
years. 

KEYWORDS: Protected areas, 
sensing, Fires, Amazonia, Climate. 


Remote 


RESUMO: Os incêndios promovem grandes 
distúrbios às florestas na Amazônia. A sua 
causa está relacionada à interação entre 
o clima e as mudanças de uso e cobertura 
do solo. Uma ferramenta muito útil nesse 
processo é a quantificação de dados 
espaciais, que através do mapeamento 
pode fornecer informações para o 
planejamento de políticas de prevenção. 
Para isso, utilizou-se dados de uso e 
cobertura do solo e desmatamento foram 
obtidos do projeto MapBiomas para a APA 
Triunfo do Xingu entre 1985 a 2019. As 
variáveis meteorológicas e focos de calor 
foram obtidas do Instituto Nacional de 
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Meteorologia e BDQueimadas. Com os resultados, a perda de cobertura vegetal durante 30 
anos foi de aproximadamente 4030,73 km2 (23,96% do território) e para 2004-2007 de 6,21% 
(1044,84 km2), 2008-2011 de 3,62% (608,52 km2), 2012-2015 de 1,40% (236,24 km2). 
Assim podemos concluir que os incêndios estão diretamente relacionados a substituição 
da cobertura vegetal por pastagens na APA Triunfo do Xingu gerando grande impacto nas 
perdas de biomassa. As mudanças de uso do solo também foram responsáveis por uma 
alteração na temperatura do ar, aumentando no decorrer dos últimos anos. 
PALAVRAS-CHAVE: Áreas protegidas, Sensoriamento Remoto, Focos de calor, Amazônia, 
Clima. 


INTRODUCTION 


Changes in land use of patterns are some of the factors that cause changes in climate 
in tropical regions (WARD et al., 2014; SPRACKLEN; GARCIA-CARRERAS, 2015; ZEMP et 
al., 2017). Deforestation and fires, for example, are present issues that directly impact the 
tropical climate (AYALA et al., 2016), such as increased temperature and reduced rainfall, 
causing greater sensitivity in ecosystems (SEDDON et al., 2016; GRIFFITHS et al., 2018). 

In Amazon, the impact of forest loss related to deforestation can be of a similar 
magnitude to the warming effects of greenhouse gases (PIELKE et al., 2016). With the 
absence of forest cover, the processes that occur on the terrestrial surface are altered, 
including lower rates of evapotranspiration and the regulation of albedo on the ground, 
which affect the magnitude and form of energy transference to the atmosphere (LEJEUNE 
et al. ., 2014; DEVARAJU et al., 2015; DEBORTOLI et al., 2016; SILVA OLIVEIRA et al., 
2018). 

The proportion of latent and sensible heat is balanced by extensive forest areas, 
which are essential in maintaining the hydrological cycle, cloud formation and precipitation 
(SYKTUS; MCALPINE, 2016). Land use changes in the Amazon have become a global 
concern, mainly due to the frequency and intensity of forest fires affecting thousands of 
hectares of forests, being responsible for the emission of large amounts of carbon into the 
atmosphere and for losses in the stock of biomass. (SALES et al. 2015; SENIOR et al., 
2017; SCOOT et al., 2018; CAMPOS et al. 2020). 

Fires promote negative health impacts, economic losses and a threat to biodiversity 
(ANDRADE FILHO et al. 2017). An increase in the occurrence of fires in the Amazon 
would further worsen the scenario and compromise policies that seek to reduce the human 
contribution to climate change (SILVESTRINI et al. 2011). This is due to the fact that fire 
is widely used in the maintenance of pastures and in agricultural areas, which when it gets 
out of control, advance into the underlying forests, reaching large areas in years of extreme 
drought (TASKER; ARIMA, 2016; ARAGÃO et al. . 2018). 

Fire usually reaches the forest from fires carried out to convert the forest into 
agriculture or pasture and to control weeds, but almost always outside the borders of interest 
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(LAWRENCE et al., 2015; SCHMIDT; ELOY, 2020). Furthermore, with the construction 
of new roads, remote areas of forest are accessed to expand slash-and-burn agriculture, 
forestry and livestock farming, consequently causing an increase in the potential conditions 
for fires to occur (SCOTT et al., 2018). 

In view of this perspective, this work analyzes landscape changes, forest fires and 
their influence on meteorological conditions observed in the Triunfo do Xingu Environmental 
Protection Area between 2005 and 2015, considering the factors of anthropic and natural 


origin in temperature variation. 


MAIN OBJECTIVES 


The objectives of this study were (i) to map the landscape dynamics in the Triunfo 
do Xingu Environmental Protection Area (APA) for 30 years; (ii) Monitor the deforestation 
behavior and air temperature variation between 2005 and 2015. 


MATERIALS AND METHODS 


Study area 


The APA Triunfo do Xingu is located in the State of Pará, created through State 
Decree in 2.612, in 4 of december of 2006, with a total area of 1,679,280.52 hectares with 
65% of its area in the territorial portion of the municipality of São Félix do Xingu and 35% 
in the municipality of Altamira (ITERPA, 2017). It faces a great management challenge, 
among other factors, because it was created in an area of high population concentration, 
with great anthropic pressure (Figure 1), with a very degraded and altered territory, and 
many land conflicts, compared to the others. conservation units created in the Terra do Meio 
region (COSTA, 2013; IDEFLORBIO, 2017). 
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Figure 1: Study area location with percentage of vegetation cover based on the Global Forest Change 
Dataset (HANSEN et al., 2013). 


Dynamics of Land Use and Coverage 


An important and promising dataset for mapping agricultural and pasture areas 
and, consequently, vegetation cover is the MapBiomas project (MAPBIOMAS, 2020). 
MapBiomas provides Landsat-based maps of land use and land cover aggregated into 5 
broad categories (and up to 25 detailed categories) from 1985 to 2019. In APA Triunfo do 
Xingu, classes are divided into: forest formation, grassland formation, pasture and water 
bodies. 

From the land use and cover change database, a historical mapping of changes 
was carried out in the period from 1985 to 2015, and consequently, the creation of a more 
detailed map characterizing the areas of deforestation that encompass different periods: 
2004 to 2007, 2008 to 2011, 2012 to 2015, and until 2019. Thus, enabling the identification 
and quantification of forest loss in each analysis period.. 


Hot spots 


The data on hotspots were obtained from the Center for Weather Forecasting and 
Climatic Studies and the National Institute for Space Research (CPTEC, 2016; INPE, 
2018) by the product called Burn Data Bank (BDqueimadas). Currently, CPTEC uses 31 
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environmental satellites (polar orbit and geostationary) to compose its observation network 
in South America (AS). The series of weather satellites are: NOAA, GOES, AQUA (EOS 
PM-1), TERRA (EOS AM-1), METEOSAT, ATSR and TRMM. These satellites perform 
orbital imaging across the country in the morning, afternoon, night and dawn. The time 
series of hotspots obtained from the Amazon corresponds to the period from 2005 to 2015. 


RESULTS AND DISCUSSION 


Change of use and coverage in APA Triunfo do Xingu 


Figure 2 provides a summary of land use and land cover changes in the Triunfo 
do Xingu APA represented by four-year highlighted images (1985, 1995, 2005 and 2015). 
It is observed that the change in land cover in APA Triunfo do Xinigu from 1985 to 2015 
was extensive. A total of 4030.73 km2 (23.96% of the territory) has undergone intense 
changes in the last 30 years. MapBiomas land use and land cover maps demonstrated the 
growing dynamics of human activities in the region that occurred from 2005 onwards, whose 
landscape changes included the expansion of pastures over the entire territory, as well as 
deforestation and forest fragmentation in protected areas. adjacent areas and indigenous 
lands (Figure 3). 
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Figure 2: Map of the dynamics of use and coverage in APA Triunfo do Xingu for (a) 1985, (b) 1995, (c) 
2005, and (d) 2015. 


Considering the period between 2005 and 2015, the territory of the APA of Triunfo do 
Xingu showed a significant reduction in the percentage of the growing number of pastures 
with the implementation of the PPCDAm (Figure 3), demonstrating the effectiveness of 
inspection to address the anthropic causes of deforestation long-term. Among the areas of 
loss, the percentage of loss of vegetation cover in the year 2004 to 2007 was 6.21% (1044.84 
km2), in 2008 to 2011 it was 3.62% (608.52 km2), in 2012 to 2015 was 1.40% (236.24 km2), 
meaning a gradual reduction of native forests between 2005 and 2015 (Figure 4). However, 
with the end of PPCDAmM in mid-2015, the loss of vegetation cover intensified until 2019 
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with a loss of 6.32% (1063.72 km2), showing the importance of inspection with articulated 


actions for territorial land planning, environmental monitoring and control and promotion of 
sustainable production activities (WEST; FEARNISDE, 2021). 
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Figure 3: Difference between periods of loss of vegetation cover and gain of pasture in APA Triunfo do 
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Figure 4: Characterization of deforested areas during the PPCDAm (2005 — 2015). 
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Changes in air temperature patterns and quantification of hotspots 


tis observed that in Figure 5, the orbital monitoring of hotspots recorded a higher 
occurrence between the years 2005, 2007 and 2015. The temporal variability of hotspots 
ends up not being a continuous phenomenon, with years with greater concentration than in 
others, and the APA do Triunfo do Xingu, described by COSTA et al. (2017) is a territorial 
extension very pressured by human actions of all kinds, such as logging, extensive livestock, 
among others. 

The year with the highest number of hotspots was also characterized by the greatest 
variation in air temperature, with 28.4ºC, compared to the climatological average of 27.3ºC. 
One of the local factors that may have influenced this increase in temperature is the lack 
of vegetation cover, which influences the prevailing weather conditions through the direct 
absorption and reflection of incident solar radiation, proven in several studies in the tropics 
(LORENZ; PITMAN, 2014; SWANN et al., 2015; HANIF et al., 2016; PITMAN; LORENZ, 
2016; LOVEJOY; NOBRE, 2018). 
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—O— Observado (ºC) 
1981-2010 (ºC) 
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2005 2007 2009 2011 2013 2015 
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Figure 5: Occurrences of hotspots and air temperature variation, APA Triunfo do Xingu between 2005 
and 2015. 


Itis also noticeable that one of the most important governance instruments at the 
federal level in Brazil, the Action Plan for the Prevention and Control of Deforestation in the 
Legal Amazon (PPCDAm) created in 2004, helped to promote a reduction in deforestation 
(Figure 4).), playing a significant role in containing environmental degradation in conservation 
units (MMA, 2013; ARIMA et al., 2014; ASSUNÇÃO et al., 2015; CUNHA et al., 2016). 
Authors such as Santoso et al. (2017) and Jiménez-Mufoz et al. (2016) also highlight 
that human actions associated with precipitation mechanisms influence the local climate, 
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as observed in 2015, during the occurrence of the El NiÃo event, classified as strong. This 
event was capable of causing scarcity of rain, low air humidity and high temperatures, 
since such action, associated with an increasing rate of hotspots and intense vegetation 
suppression, contributed to an increase in the annual variation of the temperature above 
normal, with 28.3ºC over APA Triunfo do Xingu. 


CONCLUSION 


In this study, it was observed that the association of anthropogenic factors contributed 
to a high annual variation of air temperature in APA Triunfo do Xingu. This demonstrates 
that the influence of human actions causes an acceleration of impacts at an increasing rate 
with the removal of vegetation cover if preventive actions and monitoring of deforestation 
are not implemented, such as the PPCDAm. 

Deforestation probably led to a change in the hydrological cycle of the protected 
area, causing a strong reduction in rainfall, explaining the high temperatures in 2005 and 
2007 that were not influenced by precipitation mechanisms. 

Therefore, the recycling of moisture must be considered as a key ecosystem service 
of the Amazon rainforest, in the maintenance of the local climate. 
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RESUMO: A Caatinga é considerada 
mundialmente como um ecossistema 
singular, mesmo diante de características 
climáticas adversas como o regime 
pluvial anual irregular. Aliado à esta 
situação, tal bioma vem sendo degradado 
progressivamente pela ação antrópica 
desde o período colonial do País sem 
práticas adequadas ao desenvolvimento 
sustentável. Nesse âmbito, o objetivo 
desse estudo foi apresentar uma forma 
quantitativa de avaliação do grau de 
degradação da cobertura vegetal existente 
em uma área de Caatinga em Pernambuco, 
especificadamente, na microbacia do 
Rio Terra Nova (de aproximadamente 
562.647 ha) para o ano de 2019; a partir 
de uma correlação linear entre o albedo 
de superfície, o NDVI e a Tasseled Cap 
Transformation (TCT), todos originados de 
dados do satélite Sentinel-2A. Por meio da 
análise dos resultados, pôde-se notar que 
houve a ocorrência de duas correlações 
negativas expressivas: uma entre o albedo e 
o Wetness (R? = 0.6982) e outra relacionada 
ao Brightness e o Wetness (R? = 0.7317), 
o que reforçou a influência da precipitação 
sobre a vegetação da Caatinga para com o 
processo de degradação ambiental. Para 
ambas as correlações, foram gerados os 
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seus respectivos índices de degradação de cobertura vegetal, os quais foram capazes de 
delimitar a área de estudo em regiões com degradação inexistente (21,43% e 21,02%); 
degradação baixa (36,98% e 37,78%), degradação média (28,64% e 28,37%) e degradação 
alta (12,13% e 12,01%), respectivamente. Quanto à acurácia estatística, tais mapeamentos 
alcançaram consecutivamente os índices Kappa de 0,8427 e 0,8496 e as exatidões 
globais de 88,20% e 88,72%, o que demonstrou que ambos os índices apresentam uma 
boa aplicabilidade para estudos futuros ligados à análise do nível de degradação ambiental 
existente em áreas de Caatinga. 

PALAVRAS-CHAVE: Caatinga; Tasseled Cap Transformation; Albedo de superfície. 


QUANTITATIVE EVALUATION OF VEGETABLE COVERAGE DEGRADATION 
THROUGH BIOPHYSICAL PARAMETERS FOR PERNAMBUCANO SEMIARID 


ABSTRACT: Caatinga is considered worldwide as a unique ecosystem, even in the face of 
adverse climatic characteristics such as irregular annual rainfall. Allied to this situation, this 
biome has been progressively degraded by anthropic action since the country's colonial period 
without adequate practices for sustainable development. In this context, the objective of this 
study was to present a quantitative way of evaluating the degree of degradation of the existing 
vegetation cover in an area of Caatinga in Pernambuco, specifically, in the Terra Nova River 
microbasin (approximately 562,647 ha) for the year 2019; from a linear correlation between 
surface albedo, NDVI and Tasseled Cap Transformation (TCT), all sourced from Sentinel- 
2A satellite data. Through the analysis of the results, it could be noted that there were two 
significant negative correlations: one between albedo and Wetness (R? = 0.6982) and another 
related to Brightness and Wetness (R? = 0.7317), which reinforced the influence of precipitation 
on the vegetation of the Caatinga towards the process of environmental degradation. For 
both correlations, their respective vegetation cover degradation indexes were generated, 
which were able to delimit the study area in regions with inexistent degradation (21.43% and 
21.02%); low degradation (36.98% and 37.78%), medium degradation (28.64% and 28.37%) 
and high degradation (12.13% and 12.01%), respectively. As for statistical accuracy, such 
mappings consecutively reached Kappa indices of 0.8427 and 0.8496 and global accuracy 
of 88.20% and 88.72%, which showed that both indices have a good applicability for future 
related studies to the analysis of the level of environmental degradation existing in areas of 
Caatinga. 

KEYWORDS: Caatinga; Tasseled Cap Transformation; Surface albedo. 


11 INTRODUÇÃO 


A Caatinga é o tipo de vegetação predominante na região semiárida nordestina do 
Brasil e, caracteriza-se por apresentar uma ampla variedade de espécimes florísticos e 
faunísticos unicamente encontrados nesse bioma, mesmo diante de uma irregularidade 
pluviométrica anual acentuada (SILVA et al, 2020). Em contrapartida, tal ecossistema é 
ameaçado historicamente pela degradação ambiental causada, principalmente, pela ação 
antrópica, através de práticas como o extrativismo para subsistência ou a expansão das 
áreas agrícolas (ALBUQUERQUE, LOMBARDI NETO & SRINIVASAN, 2001; RIBEIRO 
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et al., 2015; LUCENA, WANDERLEY & NÓBREGA, 2016). Fato este, que está gerando 
efeitos nocivos ao meio ambiente cada vez mais intensos devido às mudanças climáticas 
ocorrida em todo o mundo (RIBEIRO et al., 2016). 

Entre tais processos de degradação ambiental, o mais comumente encontrado em 
regiões áridas e semiáridas como a Caatinga é a desertificação, que pode ser caracterizada 
como o conjunto de alterações ambientais (diminuição e/ou alteração da vegetação nativa, 
aumento da erosão do solo, salinização do solo, etc.) que propiciam a escassez dos 
recursos naturais presentes em um determinado local (SOUZA, MENEZES & ARTIGAS, 
2015; MARTINS et al., 2019). Contudo, diversos são os exemplos encontrados na literatura 
de estudos que correlacionam a eficácia da discriminação do processo de desertificação 
em uma escala espaço-temporal através da análise de dados oriundos do sensoriamento 
remoto (CHAVES et al., 2015; LAMCHIN et al., 2016; SILVA et al., 2018; SILVA et al., 
2020). 

Em específico, o emprego de índices de vegetação aliados ao albedo de superfície 
como indicador do grau de desertificação de uma determinada área tem demonstrado ser 
uma eficaz ferramenta de análise deste tipo de degradação ambiental em regiões áridas 
e semiáridas (MA et al., 2011; PAN; LI, 2013; LAMQADEM, SABER & PRADHAN, 2018; 
WEI et al., 2018; SILVA et al., 2020). Visto que, a diminuição da cobertura vegetal acarreta 
na maior exposição do solo à radiação solar e consequentemente no aumento da sua 
capacidade de reflexão dos raios solares, expressando assim, uma correlação negativa 
entre esses dois tipos de alvos superficiais (MA et al., 2011; LAMQADEM, SABER & 
PRADHAN, 2018). 

Nesse âmbito, Lamgadem, Saber & Pradhan (2018) propuseram um modelo de 
índice de desertificação baseado na relação linear existente entre o albedo de superfície, o 
Tasseled Cap Transformation (TCT) e o NDVI, o qual demonstrou-se ser bastante eficiente 
para discriminação do processo de degradação em áreas de clima árido e semiárido. 

Portanto, o objetivo desse estudo foi avaliar a eficácia dessa metodologia para a 
determinação quantitativa do grau de degradação de cobertura vegetal ocorrido em uma 
área composta pelo Bioma Caatinga para o ano de 2019, visando à utilização desta em 
estudos futuros ligados a essa temática ambiental. Visto que, tal vegetação também se 
encontra inserida no clima semiárido, mesmo apresentando características edafoclimáticas 
diferentes em comparação a outras regiões de clima árido e semiárido ao redor do mundo. 


21 MATERIAL E MÉTODOS 


A região analisada no estudo abrange uma área de aproximadamente 562.647 ha 
e está localizada entre as coordenadas 7º45'00” S e 8º30'00” S a 38º45'00” O e 39º30'00” 
O, sendo assim, pertencente à microbacia hidrográfica do Rio Terra Nova - Unidade de 
Planejamento Hídrico, de acordo com a Agência Nacional de Águas (ANA, 2017) (Figura 
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1). Essa bacia é composta principalmente por três municípios: Terra Nova, Salgueiro e 
Cabrobó. Esta região apresenta pluviosidade e temperatura médias anuais variantes de 
400 a 750 mm e 24 a 26º C, respectivamente (ALVARES et al., 2013; FUNK et al., 2015). 


39º30W 39º15W 39º0W 38º45W 


39º30W 39º15W 39º0W 38º45W 
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' Elaborado por: Autores, 2021 
EM Brasi Fonte: ANA (2017), IBGE (2020) e Natural Earth (2021) 


FIGURA 1. Localização espacial da área de estudo. 


Os dados orbitais foram extraídos da plataforma digital do United States Geology 
Survey — USGS (https://earthexplorer.usgs.gov/), onde foram filtradas quatro imagens 
originadas do satélite Sentinel-2A referentes à localização da área de estudo (cenas 
T24MWS, T24LWR, T24MVS e T24LVR) entre os dias 05 e 15 de novembro de 2019. 
Os requisitos para a escolha desses dados foram o percentual de cobertura de nuvens 
e período de escassez hídrica (estágio temporal em que a região estaria mais propensa 
a mudanças ambientais, como as alterações da cobertura vegetal e das características 
físico-hídricas do solo), como proposto por Lamqgadem, Saber & Pradhan (2018) e Silva et 
al. (2020). 

Todo o processamento digital dos dados ocorreu nos softwares QGIS (versão 3.10.9), 
SAGA (versão 2.3.2) e Microsoft Excel 2010, e consistiu das etapas de pré-processamento 
digital, cálculo dos parâmetros biofísicos analisados, obtenção do índice quantitativo de 
degradação e análise da acurácia dos resultados encontrados no estudo. 

Primeiramente, todas as imagens orbitais passaram pelo procedimento de 
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correção radiométrica e atmosférica no plugin SCP (Semi-Automatic Classification Plugin), 
complemento específico do software QGIS, para conversão dos valores brutos desses 
dados orbitais (oriundos em números digitais - ND) em valores de refletância de superfície 
pelo método DOS — Dark Object Substraction (CHAVEZ JR, 1996; CONGEDO, 2020). 

Posteriormente, as bandas B11 e B12 foram redimensionadas espacialmente para 
a dimensão de 10 m, visando à consistência dos dados para a análise dos resultados 
(LAMQADEM, SABER & PRADHAN, 2018; SHI; XU, 2019). 

Foram utilizados os seguintes parâmetros biofísicos nesta análise geoespacial: 
Albedo de superfície (eg. 1) e NDVI (eg. 2). A equação para o albedo de superfície consistiu 
em uma adaptação do método de Liang (2001) baseada em estudos recentes (NAEGELI 
et al., 2017; LAMQADEM, SABER & PRADHAN, 2018). O NDYVI foi proposto por Rouse et 
al. (1973) e caracteriza-se principalmente por apresentar alta sensibilidade a mudanças 
sazonais e interanuais no desenvolvimento e na atividade da vegetação, como também, a 
variações do substrato presente sob o dossel foliar (JENSEN, 2011). 

Aliando essa capacidade do NDVI ao fato do albedo ser um importante indicador 
para a determinação do balanço de energia e das mudanças microclimatológicas, tem-se 
aí um meio eficaz de detectar o grau de degradação da cobertura vegetal ocorrido em 
áreas áridas e semiáridas (MA et al., 2011; LAMCHIN et al., 2016; LAMQADEM, SABER & 
PRADHAN, 2018). 


— (0.3569,+0.1300,+0.373p8+0.085p11+0.072012—0.018) 


Lsup = 1.016 (1) 
NDVI| = £8-24 (2) 
Petpa 
Em que: 


Ep — Albedo de superfície; 

Po Pp Py Pe P, Refletâncias das bandas espectrais B2, B4, B8, B11 e B12 do 
satélite Sentinel-2A, respectivamente. 

A Tasseled Cap Transformation (TCT) consiste em um índice de vegetação global 
que foi desenvolvido por Kauth e Thomas (1976) a partir de imagens do Landsat 5/MSS, 
com o intuito de produzir uma transformação ortogonal dos valores dos pixels de uma 
imagem orbital e subdividi-los em quatro espaços dimensionais: índice de brilho do solo 
(Brightness), índice do vigor da vegetação (Greenness), índice do material amarelo (Yellow 
stuff), e inespecífico (None-such). 

De acordo com o que foi proposto Lamqgadem, Saber & Pradhan (2018) e Shi & 
Xu (2019), foi gerado um TCT baseado em três indicadores espectrais: Brightness (índice 
do brilho do solo e ausência da vegetação), Greenness (índice do vigor e de diferentes 
coberturas da vegetação) e Wetness (índice da umidade do solo). Os coeficientes 
necessários a essa transformação linear encontram-se expressos na Tabela 2 (SHI; XU, 
2019). 


Ciências exatas e da terra: teorias e princípios Capítulo 5 


52 


Coeficientes do Tasseled Cap Transformation (TCT) 


TCT P, P; P, P; Pn Pa 
Brightness — 0.3510 0.3813 0.3437 0.7196 0.2396 0.1949 
Greenness — -0.3599 -0.3533 -0.4734 0.6633 0.0087 -0.2856 

Wetness 0.2578 0.2305 0.0883 0.1071 -0.7611 -0.5308 


Em que: O, D;, DP, Py P, € Pp” refletâncias das bandas espectrais B2, B3, B4, B8, B11 e B12 do 


satélite Sentinel-2A, respectivamente. 


TABELA 1. Coeficientes para o cálculo do Tasseled Cap Transformation (TCT) 


Como proposto por Chen et al. (2013) e Lamgadem, Saber & Pradhan (2018), 
para que houvesse a homogeneização dos dados, todas as imagens geradas — Albedo 


de superfície (c.,,,), 


NDVI (N), Brightness (B), Greenness (G) e Wetness (W) — passaram 


por um processo algébrico de normalização onde foi priorizada a minimização da variância 
entre tais valores, de acordo com as equações (3), (4), (5), (6), e (7): 


q = Ent 100% (3) 


A máx) Amin) 


NDVI-NDV! min) 


= E TN mm) 0 
N NDVI máx) NDVIçmin) is 100 ho (4) 
B-Bemí 
B=—— EM 100% (5) 
B(máx)- B(miín) 
G-Gemí 
=—— eim 100% (6) 
G(mãáx)—G(miín) 
W>W min) 
W = — mm 100% 7 
Wma) W min) : (7) 


Baseando-se nas definições expostas por Lamqadem, Saber & Pradhan (2018) 
em seus estudos, posteriormente a essa etapa, foram delimitados 5577 pontos regulares 
espaçados entre si por uma distância de 1000 m, os quais foram utilizados como dados 
amostrais para a análise de correlação linear (y=ax+b) existente entre os parâmetros 


biofísicos estudados. Foram determinados como variáveis independentes o Albedo e o 
Brightness, e como variáveis dependentes, o NDVI, o Greenness e o Wetness. 


Tal análise foi desempenhada através do software Microsoft Excel 2010, em que 
se objetivou investigar o nível de interação entre tais parâmetros, especificadamente em 
relação ao grau de degradação de cobertura vegetal e para tanto foi usado como indicador 
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estatístico o coeficiente de determinação entre valores ajustado (R?), como proposto por 
Ma et al. (2011) e Lamqadem, Saber & Pradhan (2018). 

Diante do que foi exposto por Verstraete & Pinty (1996), Ma et al. (2011) e Lamqadem, 
Saber & Pradhan (2018) em seus estudos, é um fato que há uma interação inversamente 
proporcional entre o albedo de uma superfície e o índice de vegetação quando se está 
analisando o grau de degradação da cobertura vegetal de um determinado local. 

Apontando assim, de forma análoga, uma correlação linear que pode ser resumida 
à seguinte expressão: 


ID=KxIV-Gap (8) 


Onde: 

ID — índice de degradação da cobertura vegetal do local; 
K— Constante encontrada a partir regressão linear; 

IV— índice espectral de vegetação; 

Gp — Albedo de superfície. 

Assim, os parâmetros biofísicos que apresentaram as maiores correlações negativas 
de forma linear foram submetidos a uma análise de regressão linear como proposto 
por Ma et al. (2011) e Lamqadem, Saber & Pradhan (2018), almejando a obtenção da 
equação pertinente ao grau de degradação existente na área de estudo para o ano de 
2019. Posteriormente, foram gerados os mapas de índice de degradação a partir dessas 
equações criadas, obedecendo a quatro níveis: Degradação inexistente, degradação baixa, 
degradação média, e degradação alta (MA et al., 2011; LAMQADEM, SABER & PRADHAN, 
2018; WEI et al., 2018). 

Para a segmentação das classes foi utilizado o método de quebras naturais de 
Jenks, que consiste no agrupamento natural inerente aos dados, sendo seu limite definido 
pela posição onde os valores dos dados tornam-se relativamente diferentes e foi escolhido 
por ter apresentado alta acurácia em diversos estudos lidados a análise do grau de 
desertificação (MA et al., 2011; PAN; LI, 2013; LAMQADEM, SABER & PRADHAN, 2018; 
WEl et al., 2018). Tal método é capaz de maximizar as diferenças entre as classes geradas, 
agregando efetivamente os valores com mínima variância entre si em classes semelhantes, 
de modo a obter um ótimo resultado de classificação (ZENG-HAI; TAO, 2013; WEI et al., 
2018). Foram utilizados os mesmos 5577 pontos amostrais para a execução de tal etapa. 

Visando-se impedir a interferência nos resultados, foi necessário aplicar uma 
camada de máscara para a remoção das áreas urbanas e dos corpos hídricos pertencentes 
à região analisada nesse estudo, como proposto por Lamgadem, Saber & Pradhan (2018) 
e Wei et al. (2018). 

Para a análise da acurácia, primeiramente foi criado um mapa de uso e cobertura 
de solo por meio de uma classificação não supervisionada obtida pelo uso algoritmo 
K-means (lotado no software SAGA — versão 2.3.2), de dados do Projeto MAPBIOMAS e de 
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interpretação visual de imagens de satélite, onde foram extraídas e delimitadas as classes 
de cobertura vegetal para a região de estudo que serviram como referência para validação 
do índice de degradação produzido (LAMQADEM, SABER & PRADHAN, 2018; WEI et al., 
2018; MAPBIOMAS, 2021). 

Posteriormente, foi escolhido um total de 390 pontos amostrais aleatórios para a 
realização do cálculo da acurácia dos resultados encontrados a partir de matriz de confusão, 
especificadamente, a partir dos seguintes parâmetros estatísticos: índice Kappa e exatidão 
global (MA et al., 2011; LAMQADEM, SABER & PRADHAN, 2018). 


31 RESULTADOS E DISCUSSÃO 


A partir da análise de correlação linear realizada no conjunto de dados normalizados 
utilizados nesse estudo, notou-se que a maior interação entre os parâmetros foi encontrada 
entre os índices Brightness e Wetness (R? = 0.7317), seguida pela correlação entre o albedo 
e o índice Wetness, com um coeficiente de determinação de 0.6982 (Figura 2A e 2B). Por 
outro lado, os índices NDVI e Greenness não apresentaram significância na correlação, 
ambos com valores inferiores a 0,30. Resultados semelhantes foram detectados por 
Lamqadem, Saber & Pradhan (2018) em seus estudos, que enfatizam a relação existente 
entre a diminuição da umidade do solo e o processo de desertificação ocorrido em tais 
áreas áridas e semiáridas: a redução da água no solo propicia a aceleração do processo 
de degradação ambiental através de fatores como a alteração da cobertura vegetal, 
aceleração do intemperismo, o acúmulo de sais minerais no solo (SANTOS et al., 2017; 
MARTINS et al., 2019). Isto torna essas áreas mais propensas ao processo de degradação 
em si (LUCENA, WANDERLEY & NÓBREGA, 2016; SILVA et al., 2018; SILVA et al., 2020). 


y=-0,9143x + 85,936 y=-0,9071x + 86,066 
100 o é Rº=0,7317 100 o Rº2=0,6982 
Ba 99 O da O 


Wetness 
Wetness 


0 20 40 60 80 0 20 40 60 80 
(A) Brightness (B) Albedo 


FIGURA 2. Análise de correlação linear: Brightness-Wetness (A) e Albedo-Wetness (B). 
Assim, as constantes K encontradas foram 1,0937 e 1,1024 (Brightness-Wetness 
e Albedo-Wetness, consecutivamente), e através da técnica de quebras naturais (Jenks), 


pôde-se realizar a classificação da área de estudo de acordo com o nível de degradação 
(inexistente, baixo, médio e alto), a qual está representada nas Tabelas 2 e 3, tanto para a 
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correlação Albedo-Wetness quanto para Brightness-Wetness, respectivamente. 


Área Área índice Exatidão 
Classe de degradação ID mapeada mapeada Kabpa Global 
(ha) (6) Er (%) 
Degradação inexistente ID > 52,23 120.565,15 21,43 
Degradação baixa 38,46 <ID<52,23 208.064,66 36,98 
ss 0,8427 88,20 
Degradação média 20,78< ID < 38,46 161.166,59 28,64 
Degradação alta 20,78 >ID 68.254,84 12,13 


TABELA 2. Classificação de acordo com o índice de degradação da cobertura vegetal para a correlação 
Albedo-Wetness. 


Área Área índice Exatidão 
Classe de degradação ID mapeada ' mapeada Kappa Global 
(ha) (6) ii (%) 
Degradação inexistente ID > 52,3 118.284,2 21,02 
Degradação baixa 38,07<ID<52,3 212.562,06 37,78 
; 0,8496 88,72 
Degradação média 20,96 < ID < 38,07 159.638,26 28,37 
Degradação alta 20,36 > ID 67.566,44 12,01 


TABELAS. Classificação de acordo com o índice de degradação da cobertura vegetal para a correlação 
Brightness-Wetness. 


Nessas tabelas, pôde-se notar que em ambos os casos as classes que se 
destacaram foram as de degradação baixa e média, com um total de aproximadamente 
208.064,66 e 161.166,59 ha, respectivamente, para a correlação Albedo-Wetness e de 
212.562,06 e 159.638,26 ha para a correlação Brightness-Wetness. Contudo, vale ressaltar 
que, mesmo sendo estas as maiores classes entre as demais, de modo geral, ainda houve 
um alto quantitativo de áreas que estavam sujeitas ao nível alto de degradação, totalizando 
68.254,84 ha (correlação Albedo-Wetness) e 67.566,44 ha (correlação Brightness- 
Wetness). Quanto à acurácia dos mapeamentos realizados, foram encontrados os índices 
Kappa de 0,8427 e 0,8496 e as exatidões globais de 88,20% e 88,72%, consecutivamente; 
demonstrando assim, a ótima eficácia de detecção da metodologia proposta para avaliar o 
processo de degradação da cobertura vegetal ocorrido na área de estudo no ano de 2019 
(MA et al., 2011; LAMQADEM, SABER & PRADHAN, 2018). 

A Figura 3 apresenta uma discriminação espacial quanto aos diferentes níveis de 
degradação de cobertura vegetal existentes na microbacia do Rio Terra Nova, segundo os 
respectivos índices quantitativos de degradação gerados neste trabalho. 
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FIGURA 8. Classificação quanto ao nível de degradação da cobertura vegetal para a área de estudo 
durante o ano de 2019. 


41 CONSIDERAÇÕES FINAIS 


A proposta desse estudo foi apresentar uma forma de avaliação quantitativa que 
fosse capaz de detectar os diferentes níveis de degradação de cobertura vegetal existentes 
em áreas da região semiárida do Nordeste Brasileiro. Assim, diante dos resultados 
encontrados após a validação dos dados, pôde-se concluir a eficácia desta metodologia para 
a determinação do processo de degradação da cobertura vegetal ocorrido espacialmente 
em regiões de Caatinga como a área de estudo. 

Os índices gerados ainda foram capazes de detectar a correlação existente entre a 
quantidade de água presente no solo (representada pelo índice Wetness) e o percentual de 
solo exposto presente no local de análise (representado pelos índices Albedo e Brightness), 
representando assim, uma menor presença de vegetação neste local. Característica esta, 
intrínseca à Caatinga em épocas de seca. 
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